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The lowest order radiative corrections (RC) to width and spectra of radiative 7re2 
- decay are calculated. We take into account virtual photon emission contribution 

■ as well as soft and hard real photon emission one. The result turns out to be 

O 

' consistent with the standard Drell-Yan picture for the width and spectra in the 

leading logarithmical approximation which permits us to generalize it to all orders of 
perturbation theory. Explicit expressions of nonleading contributions are obtained. 
. The contribution of short distance is found to be in agreement with Standard Model 

> ; 

1/-^ ' predictions. It is presented as a general normalization factor. We check the validity 

CSl ' of the Kinoshita-Lee-Nauenberg theorem about cancellation in the total width of 

o 

the mass singularities at zero limit of electron mass. We discuss the results of the 

m ■ 

I previous papers devoted to this problem. The Dalitz plot distribution is illustrated 

Q^' numerically. 

D , 

■ I. INTRODUCTION 



^ • The process of radiative negative pion decay 



7r-(p) -^e-(r) + z7e(g)+7(fc) (1) 

attracts a lot of attention both experimentally and theoretically [H, P|, S], 1^; [3; Ql- 
The main reason is a unique possibility to extract the so-called structure-dependent part Msd 
of the matrix element 

M = MjB + Msd, (2) 

which can be described in terms of vector and axial-vector formfactors of the pion and the 
hints of a " new physics" including a possible revealing of tensor forces (which are absent in the 
traditional Standard Model). 

Due to the numerical smallness of Msd compared to the inner bremsstrahlung part Mjb, 
the problem of taking into account radiative corrections (RC) becomes essential. The lowest 



order RC for a special experimental set-up was calculated in 1991 in '4^, where, unfortunately, 
the contribution to RC from emission of an additional hard photon was not considered. This 
is the motivation of our paper. 

Our paper is organized as follows. In part m we calculate the contributions to RC from 
emission of virtual photons (at 1 loop level) and the ones arising from additional soft photon 
emission. For definiteness, we consider RC to the QED-part of the matrix element M/^ 



Mtb = -iA uJr) 



{6*r) {e*p) ^ e*k 



Gf 

[1 + '^5)vu{q), A = e-j=Vudf7rm, (3) 



(kr) [kp] 2{kr) 

where e = y/47m, Gp = 1.17 ■ 10^^ GeV~^ is the Fermi coupling constant, e is the photon 
polarization vector, f^^ = 131 MeV is the pion decay constant, Vud is the element of the 
Cabbibo-Kobayashi-Maskawa quark mixing matrix, and M and m are masses of pion and 
electron. The explicit form of the matrix element including Msd is given in Appendix 1X1 

In part lllll we consider the process of double radiative pion decay and extract the leading 
contribution proportional to "large" logarithm In M^/m^ which arises from the kinematics of 
emission of one of the hard photons collinearly to electron momentum. We arrive at the result 
which can be obtained by applying the quasireal electron method 0|. 

In conclusion, we combine the leading contributions and find that the result can be expressed 
in terms of the electron structure function j^. Really, the lowest order RC in the leading 

/ 2 \ " 

logarithmical approximation (LLA) (i.e., keeping terms ( ^ In ^ j ) turns out to reproduce 
two first order contribution of the electron structure function obeying the evolution equation 
of renormalization group (RG). This fact permits us to generalize our result to include higher 
orders of perturbation theory (PT) contributions to the electron structure function in the 
leading approximation. In conclusion, we also argue that our consideration can be generalized 
to the whole matrix element. 

As for the next-to-leading contributions, we put them in the form of the so-called i^-factor 
which collects all the non-enhanced (by large logarithm) contributions. Part of them, arising 
from virtual and real soft photon emission, is given analytically. The other part, arising from 
emission of additional hard photon in noncoUinear kinematics is presented in Appendix O in 
terms of 3-fold convergent integrals. 

In Appendix ^ we give the simplified form of RC (the lowest order ones) and make an 
estimation of the omitted terms which determine the accuracy of the simplified RC. 

Appendix O contains the list of 1-loop integrals used in calculation of 1-loop Feynman 
integrals. 
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In tables |lj [HJ II I II the result of numerical calculation of contributions of RC to Dalitz-plot 
distribution in Born level, leading and non-leading approximations are given. 

II. VIRTUAL AND SOFT REAL PHOTON EMISSION RC 

A rather detailed calculation of the lowest order RC was carried out in j^]. Nonclear ma- 
nipulations with a soft photon emission contribution used in Ql results in a wrong form of the 
dependence of RC on the photon detection threshold Ae which contradicts the general theorem 
0]. Another reason for our revision of [4J is the mixing of (on-mass shell and dimensional) reg- 
ularization schemes in it. We use here the unrenormalized theory with the ultraviolet cut-off 
parameter A.^ 

Following ^4], we first consider RC to the "largest" contribution - inner bremsstrahlung Mjb- 
As well we consider the central part of Dalitz-plot and omit (if possible) [3 = (^)^ = 1.34-10"^, 
that is /? — s> 0. 

We distinguish the contribution to RC from emission of virtual and soft additional photons 
(Ay and A^, respectively): 

E \MlBfL,^,^,, = j:msf{l + ^{Ay + As)), (4) 



spin states 



where 



^-^ x'^ix + y — 1 



x'^^x + y 

(here x = y = z = M is the pion mass). 

Soft photon emission RC have a standard form (see, for example j8|, formula (16)) 

(Pki f p - \ 2 



a . a 
-As 



jx = \/kl+X^<Ae<t.yM 



a 

TT 



(Le - 2) In + 11, - + 1 - 6 I (1 + om , (6) 



2 2 

where ki the additional soft photon momentum, = In ^ and A is the "photon mass" , ^2 = ^• 
This result agrees with the general analysis of infrared behavior given in 0]. 

Now let us consider the calculation of the virtual photon emission corrections. First, we use 
the minimal form for introduction of the electromagnetic field through the generalization of the 
derivative p^j, p^ — ieA^. 






where 





m 



FIG. 1: Effective vertex. 



The sum of contributions (all particles except antineutrino can be on or off mass shell) leads 



to: 



elf + k + m)p , . ^ , . i(f + k + m) . . 



(r + ky 



(r + kY — "rn? 



p = r + k + q. Thus, we can extract the contact photon emission vertex and obtain the effective 
vertex ~ m, as it is shown in Fig. ^ 

In terms of this new effective vertex we can write down the 1-loop Fevnman diagrams (FD) 
of virtual photon emission RC (f Ay) which are shown in Fig. |21 (see |^). We should notice 
that these diagrams can be separated out into four classes. The contributions of each of these 
classes are gauge invariant. 

Three first classes were considered in Q], where their gauge invariance and in particular 
zero contribution of class II were strictly shown. The last statement directly follows from the 
gauge-invariance . 

The contribution of class III contains the regularized mass and vertex operators. The relevant 
matrix element has the form: 



M, 



III 



-iA 



a 



27iM 



where 



Ue{r) 



1 1 



e* -k 



(e*r) 
{kr) 



[l + -i^)v^{q), 



(8) 



—n H 

y' 2a 



a 

1 ol/2 



-Is 
21/ 



3l/' + 2 



2y'a 



In2 



B2 = - 

n = l.[Lt2il)-Lt2il + y')-t7T\n{l + y') 
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Class I 




Class II 




Class III 






Class IV 



FIG. 2: Feynman diagrams of virtual contributions. 



here = In | — ivr, a = 1 + y', y' = We see that Mjjj is explicitly gauge-invariant and 
is free from infrared singularities. At the realistic limit y ^ x ^ z ^ j3 (rather far from the 
boundaries of Dalitz-plot) we obtain for contribution to the matrix element square (structure 
~ Bi gives a zero contribution in the limit (3 —>■ 0) (in agreement with 0|): 

y 



A\Mnif = 2Re{iMjB)^Mjn} = A 



a 



2 \f3 



(9) 



TT XZ 

As we work within the unrenormalized theory, we should consider class IV (not considered in 

n 

jj|) which concerns the contribution of counter-terms due to renormalization of the pion and 
electron wave functions (see formula (17)) 



A \Mjv\ 



1 3 , ^ , M2 9 
--La + -ln/3 + ln— 

2 2^ A2 4 



, M2 3' 



(10) 



where La = In and A is the ultraviolet cut-off parameter to be specified later. The first 
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term in the brackets in the r.h.s. is the electron wave function renormahzation constant and 
the second one corresponds to the pion wave function renormahzation. 

Let now consider the contributions of FDs of class I. In [J] the explicit gauge-invariance of 
the sum was demonstrated: (Tf + + + + T^) = 0, where Tj^ is the contribution of 
the corresponding FD of Class I (see in Fig. Ej). 

The contribution to the matrix element square can be written as: 



A \Mjf = 2Re {{Mjb)^ Mj] = ^A^Re y ^ [A, + A, + A, + A, + A, 



(11) 



where 



^1 = 

^3 = - Sp [r{2p - k){r - fci)(l + 75)g {Q{p - ^i)) + ~ '''' 



(0)(1)(10(2) ^ ^ ^ - '-^ y-^^ 

1 ^ ^ ^ ^ / 

A, = - (o)(i)(2)(20 ^^^^^ ~ ^'^^^ ~ ^'^^"^^ + k-h){l + 75)g I + ^ 

^5 = ~l (o)(i)(20 ^^ + ^^^^^ - fci)(r + - + 75)g (|^' + ^ j]]^- 

where /ci virtual photon momentum, = — (f^ + ~ IP' ("7 ~ Here the denominators 
of the pion and electron Green functions {i) are listed in Appendix^ We omitted m everywhere 
in the numerators. The ultraviolet divergences are present in Ai and A5. 

Using the set of vector and scalar integrals listed in Appendix O and adding the soft photon 
emission corrections we obtain 



1 + - (Ay + As) 
TT 

= <! 1 + - 



{l + ^[(L.-l)(lnA + ^)-lnA + ^..^(.,,)]} 



(12) 
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where Sw = I + |f i^A, ^ = ^- Ri^, v) looks like 

/ 3 2\ 

R{x,y) = Aahy {—blny + - — A\n.y\n-\ — 

- ^hhiy {-ay{l + ln4) + ((x - 2f - y{a + 2)) \nx + x{l + y^) In^) + 
+ 2b\n^ y {-16y + x{16 - 7x + 7y{x - 2))) - 

- 6 ■ 8by{6 + 2x^ - a;(5 + y)) + 

+ Axyzh ^Lz2 ^ — Li2(l — y)^ + 8abyLi2{x) + 

+ 8aby In xln{l — x) — 8xxyz{x + y — 2)lnz — 

- Axy z In^ z{—x + y{x — 2) + y"^) — Axyb, (13) 
B'{x,y) = 8{l~y){l + {l~xy){x + y-2)^ = 8aby, (14) 

where a = 1 + (1 — x)^, b = {x + y — 2)^, x = l — x,y = l — y, z = x + y — 1. Here we do 
not have a complete agreement with the result of Note that the sum A5 + Ay does not 
depend on "photon mass" A, and, besides, the coefficient at large logarithm agrees with RG 
predictions. 



III. EMISSION OF ADDITIONAL HARD PHOTON 

Emission of additional hard photon (not considered in the previous papers concerning RC 
to the TT — i> ez/7 decay), i.e., the process of double photon emission in the ne2 decay: 7r^(p) 
e~(r) + Ueiq) + 7(^1) + 7(^2) is described by 11 FD drawn in Fig. El 

In collinear kinematics (when one of the real photons is emitted close to the electron emission 
direction) the relevant contribution to the differential width contains large logarithms. 

The matrix element of the double radiative pion decay has the form 



j^^— 77 = ^A./^ [ne(r)0^,(l + 75)^.(g)]-£^(fci)4(fc2), (15) 

where £1,2 are the photons polarization vectors which obey the Lorentz condition: (k\ e\{k\)) = 

(^2 £2(^2)) = 0. The tensor O^'^ has the form 

^ ,u . 2p^-2{p-k,r 2p'^ ■ 2{p ~ k,Y \ 1 

V -2(M -2(M ) {p-k^-k^f-M^^ 

2pM Y{r + k2 + m) 2p'' Y{r + ki + m) 
^ -2(M 2(A;2r) ^ -2(M 2{k^r) ^ 



+ 



^f + k2 + m ,j + ki + m ^ 

^ 2{k2r) ^ ^ 2{kir) ^ 



f + ki + k2 + m 

(16) 



(r + /ci + k2y — rn? 






2 1 








1 2 




FIG. 3: Feynman diagrams, describing the double radiative pion decay. 



One can be convinced in the exphcit fulfillment of the requirements of Bose-symmetry and 
gauge-invariance M'"~"^^"'"' {ei ki) = M'^~*'^'^"'"' {e2 —>■ ^2) = 0. The expression for 
^ \M^^'^'^^^f ~ Sp [{f + m)0^^q O^y] is rather cumbersome. 
The contribution to the differential width has the form 



7-phard 



2Af 



(17) 



6^{p ~ q ~ r ~ ki — k2 



(27r)i2 2qo 2ro 2wi 2^2 

We do not take into account the identity of photons: we believe the photon with momentum ki 



to be a measurable one with 

2(fc2p) 



2(fcip) 
M2 



X and the photon with momentum ^2 to be a background 



one with -^^r^ > A and most general kinematics. 

It is convenient to consider separately the coUinear kinematics of emission of one of the 
photons, i.e., the case when the angle of emission of one of the final photon direction of motion 
is rather small 6i = (r, ki) < 60 I, 6q ^ (3. (Note that double collinear kinematics is excluded 
since the invariants x'12 = ^^^^2''^ cannot be small simultaneously). 

The contribution of these collinear kinematics to the differential width contains "large loga- 
rithms" Lf.. To extract the relevant contribution we can use the quasireal electron method 
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For this aim let us arrange the integration over phase volume in the following way: 



|M|^ c/$4 = \Mc\^ d^^ + \Mc\^ (^$4 - d^^) + (|M|^ - \Mc\^) ^$4, 



(19) 



with 



y'' + {y + x2f 10 



x'2 X2{x2 + y) 



X' 



i2 



{l + {l-xif){l-X2-y) 
x\{xi + X2 + y - I) 



(20) 



(27r)«M-d$4 = M-^'^ ^'^^ ^'^^ '-^^ 



2ro 2u;i 2cc;2 



5 ((p - r - /ci - /C2)^) 



77 



— xiX2y dxidx2dy dCi dVt2 x 



X (5 ^1 — xi — X2 — y + 



^fi(i - c„) + f^d - CO + (1 - 



(27r)«M-4rf$4^ = - 



2^ ?/ + X2 



dxidx2dy dil!. 



2 ) 



.(21) 
(22) 



where X12 



^'1,2 = C12 = cos{ki,k2), Ci,2 = cos(A;i,2,r), (ifi2 is the angular 

phase volume of the photon with 4-momentum ^2- 

Note that integration over the angular phase volume dQ2 in (see (jl9p ) is restricted by 

2n 1 

J dn2 = f d(j) f dC2 and in the second term in (fT^ (i$4 — can be replaced by c/$4 

1-^2/2 

27r l-eg/2 

with J (ii72 = J d(j) J dC2- The second and the third terms in (fT^ do not contain collinear 
-1 

singularities, i.e., are finite in the limit /? ^ 0. 

The contribution of hard photon emission can be written out in the form 

1 



J'phard 
IB 

dxdy 



t 



+ 



(23) 



+ -K^^^''{x,y,9o,A), 

IT 

where the function B{x,t) is defined in (0), 

and Kj^^'^{x, y, 60, A) represents the contribution of two last terms in ()19j] . Now it is convenient 
to introduce the following quantity: 

1 



6(1- z- A) 



(24) 



KUx,y,A) = K^^^''{x,y,9o,A) + J jB{x,t) 



, In ^ + 1 - ^ 
tJ 4 t 



(25) 
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which does not aheady depend on the parameter ^o- The exphcit view of K^q{x, y,A) is given 
in Appendix O 

In the total sum of virtual, soft, and hard photon emission contributions all the auxiliary 
parameters - "photon mass" A and A - are cancelled out. The resulting expression for the 
differential width with RC up to any order of QED PT with the leading logarithm ^L^ ~ 1 
and the next-to-leading accuracy have the form 

g^^^.„,/fB(M)I'(f)(l + ^/W...)). (26) 



where D{z) is the well-known electron structure function, which has the form 

1 / a \2 



where 



D{Z) = ^(1 - ^) + 1^ (Le - 1) + ^ (£) (Le - 1)' P^'Hz) + 



P(^) = lim 



^-t^e{l - z - A) + 6{1 - z) (2\nA + ^ 



1 + z^ 
1 - z 



(27) 



(28) 



p(")( 



^p(l)(^)p(n-l) P 



(29) 



The function KjB{x,y) is the so-called K-factor which here has the form 

Kjb{x, y) = - 2 In A + K'^^{x, y, A), (30) 

where R{x,y) and B\x,y) was defined in (fT^ and (fT^. The quantity Kj^{x,y, A) — 2 In A is 
finite at A ^ and its explicit expression is given in Appendix O 



IV. DISCUSSION 



It is easy to see that in the total decay width all the dependence on /3 disappears in accordance 
with the Kinoshita-Lee-Nauenberg (KLN) theorem j^. Really, we obtain integrating over y: 

11 It 11 1 

jdyjjD (I) ■ fit) = j dt fit) I (I) = jdt fit) I dzD iz) = I dtfit). (31) 



Now let us discuss the dependence on the ultraviolet cut-off A. It was shown in a series 
of remarkable papers by A. Sirlin that the Standard Model provides A = M^. Another 
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important moment (not considered here) is the evolution with respect to ultraviolet scale of 
virtual photon momenta from hadron scale {rrip) up to Mz It results in effective replace- 
ment Sw —>■ Sew = 1 + - In — f- ^ 1.0232. So all the QED corrections to the total width are 

n rrip 

small ~ O (^), but the electroweak ones rather large: F ^ Fq ■ Sew- The factor Sew can be 



absorbed by pion life-time constant fn\/^EW —>■ fT^ Thus, we replace A, defined in ((31), 
with Af. 



'-exp. 



A = e%VudUm ^ = e%Kd/r^m. (32) 

v2 I — I v2 

We also note here that according to jjjj we must use the redefined constant A^xp. in the Born 
approximation. That is why we use A^xp. in Appendix 1X1 

In our explicit calculations we considered RC to the inner bremsstrahlung part of the matrix 
element. Let us now argue that in the integrand of the r.h.s. of fl26j) one can replace B{x,y) 
by the total value including the structure dependent contribution $|°](x, y) which is defined in 
Appendix El in ()A5j) . This fact can be proved in the leading logarithmic approximation. One 
of contributions arising from additional hard photon emission close to the electron direction 
can be obtained by applying the quasireal electron method and has the form with 
p(i) ^ p(^\ The KLN theorem in a unique form provides the soft photon emission and virtual 
RC to be of a form with complete kernel P{z) in the structure function D{z) in (j26|l . So our 
result reads 

1 



dxdy 



I f <Hx,t)D (I) (l + ^K{x,y)) . (33) 



We also calculate the RC to the SD part in leading logarithmical approximation. Formula 
()33|) can be trusted in the region where the IB part dominates. In the region where IB < 
SD, which is suitable for SD measurement, both IB and SD are small; so the question about 
non-leading contributions becomes academical. Thus, we suggest here that K-factor is the same 
order of magnitude as for the IB part (i.e. K{x,y) = KjB{x,y)), where it can be calculated in 
the model-independent way. Its numerical value is given in Table ITTTl 

2 

We underline that the explicit dependence on " large logarithm" Le = In ^ is present in the 
Dalitz-plot distribution. 

Now let us discuss the results obtained in some previous papers devoted to RC in the 
radiative pion decay. 

In the hard photon emission was not considered which led to violation of the KLN- 
theorem. 
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In the main attention was paid to possible sources of tensor forces. As for real 

QED+EW corrections depicted in formula (13) in Q: ? ~ 0.7f ~ 0.2% the QED leading 
RC was presumably omitted. 
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APPENDIX A: BORN AMPLITUDE 

The radiative pion decay matrix element in the Born approximation has the form Q| 

M = MiB + MsD + Mr, (Al) 

where 



MiB = -iAexp. Ue{r) 



{e*r) {e*p) ^ e*k 



(kr) (kp) 2{kr) 



'l + lB)vu{q), (A2) 



MsD = fr/e.,. Ue{rh,{l + l,)vM [Py E^'^^^p^kp + iF^ {g^'' (kp) - p^kn] , (A3) 

A 



Mt = -i^^. Ft e^k" u,{r)a^,{l + ^^)vM - 

TTLJti 

-i "^llxp. Ft (^«e/3 - k/sCa) ^^" Ueir)a°"'{l + i5)vu{q), (A4) 



where Q = p — k. The constant A^xp. defined in (j32|) . Squaring amplitude ()A1|) and summing 
over final photon polarizations leads to the following decay width: 

dx dy ° (27r)' 

+2v^ [ifv + fA) ^fl^ + {fv - Ia) + 

+ 2 {2fT{fT - f^) + f?) + 2 (2(/^ - /^) + f!,x) $g) } , (A5) 
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where 



,<«..(i^|Ki±(i_£a.B,...), 

x^{x + y — 1) 
= + 1)2, 

= (l-x)(l-^/)^ (A6) 
*?nU = -zi^-x){l-y) 



X 



x^ ' \ x + y 

4°) = (l-y){x + y-l), 



(0) _ 1-1/ 



where fv,A,T,T' = {M"^ /2mfT,)Fv,A,T,T' ~ 146 ■ Fy_A,T,r'- In particular, the conservation of vector 
current hypothesis relates the vector formfactor Fy to the lifetime of the neutral pion 



1 /2r(7r0 ^77) , , , 

.^^uji = -a/ — ^ ' ' = 0.0259 ± 0.0005 (A7) 

or equivalently fv ~ 3.78. The Dalitz-plot distribution of the inner bremsstrahlung part of the 

>(0) 
IB 



Born amplitude is given in Table H] 



APPENDIX B: SIMPLIFIED FORMULA FOR RADIATIVE CORRECTIONS 



Here, we present the simplified form of radiative corrections 



RC 



+2v/^ {{fv + fA) <U + Uv - Ia) ^fL) + 
+2 {2frUT - f^) + f?) ^f, + 2 (2(/t - /^) + M < 
1) [^^il + (/v^ + fA? + (/y - fA? + 

+2v^ ((/v. + fA) $ + (/v - Ia) + 
+ 2 (2/r(/T - /^) + /^') ^^rl + 2 (2(/r - fr) + /t^^) <f g ] } , 



(Bl) 
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where ^f'^ was presented in Appendix 1X1 and 



(1) 

IB 

(1) 

+ 

(1) 



(1) 
Int+ 



(1) 
Int- 



(1) 



1 +X 



^2 



3y y x + xy , r^V . V x^x"^ + y"^) x 
-- + ^ — ln?/ + 2-ln — In- 

2 z X x^ z y x'^z z 



X 

X 

X 
X 

1 

X 



-2;^ H — h y{y — 2a;) + x(x — 1y) In y — x^y + It? In — 



\f + + - 3) + (1 - 2|/) In y + 2y2 In ^ 

y -1 0-1 

y m y — 2?/ m — 

2 y. 
1 3 x^i/ / y \ 

— -xy + h X I ?/ In ?/ + 2?/ In - ) + 



2 " 2 z 

2 



+ - - 



y x + xy , ^y ^ y x(x^ + 1/2) 
In y + 2- In 



x^ 



z y 



x'^z 



z 



^1 " 2 

(1) _ y 



-yz- ^ +y(2x + y) - (xy-y) \ny + 2yz\n-, 



X 



-- + \ny + 2ln- 
2 2/. 



(B2) 



where z = x + y — 1, x = 1— x, y = 1 —y. Here we should notice that the functions ^[^^ satisfy 
the following property: 

1 



/ 



dy $ 



(1) 



0, 



(B3) 



which is in accordance with the demands of the KLN-theorem. 

Let us now estimate the magnitude of the terms omitted in our approximate formula ()B1 
They are 



my 

m) 



a 



/m\2 a, M2 /a, M'^V 



IT 



-K{x,y) : j .-In— :(-ln— ) = 10"^ : 1 - 2 ■ 10"^ : 10"^ : 10-^(B4) 



respectively. The main error arises from which is presented in fl3()|l and Appendix O The 
Dalitz-plot distribution of RC to the inner bremsstrahlung part of the Born amplitude is 
given in Table ITTl 



APPENDIX C: HARD PHOTON EMISSION ivTj'g-FACTOR 

In numerical calculation of the A'j^^-factor ()25|) it is convenient to use the following form 
of phase volume dHJ: d^4^ = -j^wfrTlT ^y '^^2, where ^1^2 = dC2d(f)2, A2 = 2 - 
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a;(l-Ci2)-|/(l-C2),Ci,2 
emission RC reads 



cos(A;i^2, r). Thus the X^^-factor which comes from hard photon 



K^j,{x,y, A) -2ln A 



47r 



dCidQ2 



\A2 



1 



B{x,y)x'2 



h{x,X2,y) - \ A\ 



X2 + y 



+ 



dt B{x,t) 
'TB{x,y) 



(CI) 



where 



lL{x,X2,y) = -y 



_ ^ _ ^2 ^ ^ _ 111. + y^. 

a a a b 

'\ 9 9 9 

y X yx y X 

+ ^ + ^—^ +X2-yx2- 



_ yx^ 

a b 
2yx2 



b 

yx2 



a b 

y"^ X2 



ab 

2X2 



a b 



b X 

,2 



^ 2yx2 ^Ayx2 XX2 ^XX2 



2y'^x 
a b 
_ '2yx2 

a b 
yxx2 



+ 

yxx2 



X 



ax 



62 



2y^xx2 



a b 

X^ X2 



X^ X2 ^yx^ X2 ^ x-i 

&2 ab b 



+ 



XX2 



3 



62 



y^ X ^ yx 



62 
2 



x'l x'l 



yx2" 

b 

yx 



y2 x2 



yx2 

x'l 



yx 

x'l bx'i 

yx2 



+ 



yxx2 

9 

y x 



62 

2 



+ 



+ 



yxx2 



a 6 



bx'i bx'i bx'i 

2 y X X2 2 y x x'2 6 y x X2 



x'l 



2y XX2 



62 x'l 

4 y ,x2 X2 



bx'i 



bx'i 

4 y2 X X2 
62 x'l 



■u2 X2 



2 1/2 a; 
x'l bx'i 

2 1/2 X2 



y^ X2 



X'l 



bx'i bx'i 



+ 



62 x'l 

2 



4 y a; a;2 ^ Ay x X2 



6 y^ X X2 
bx'i 
2yx^X2 



yx2 



y^x2 



62 x'l bx'i 62 , x'l 62 x'l 

2 Ayxx2'^ Ayxx2^ 4:y'^xx2^ 



+ 



6 x'l 6 x'l 



62 x'l 



bx'i 



+ 



62 x'l 



+ 



Ayx^X2^ 2yxx2 



62 x'l 

2 X2 x'l 



+ 



+ 



62 x'l a ab 

2 X2 x'l 2 X2 x'l 4 X2 x'l 

ab X ax 



y"^ xx'i ^2x2 x'l ^ 



a 6 



a 



XX2X1 XX2X1 



62 



a 6 



2yxx2x'i xx2^x'i 



a 6 



a 6 



(C2) 



W\ T singular / \ , rregular / \ 

x,X2,y) = — /jvc {X:X2,y) + lNc [x,X2,y), 



X2 



(C3) 
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jsingular / \ A X Ay X Ix lyX AljX Ax X i 

^NC l^;, X2,2/j — I 5 - I - I - 5 

a x \ x \ ax I 

Axx'o 2xx''2 Axx'o 



or x'l ax' I 

jreguiar, \_o 8 8 y _ A y Q X A X _ Ayx 3yx _ X _ x'^ 

a a a ab ao ao ab 

yx^ 6x2 4x2 4?/X2 3yx2 y^ X2 Ax2 Ayx2 



(C4) 



+ + ^ + + — + 



ab a ab a^ ab ab ax a^ x 

2yxx2 6yxx2 3x'^X2 X2 yx2 'ixx2 y ^ 

a^ ab ab ab ab ab x'\ 

1y y"^ 2y'^ ^ X yx 2yx yx 2yx ^ 



a x'l x'l a x'l x'l x'l ax'i bx'i abx'i 
y"^ X yx"^ 2yx2 yx2 yx2 y^ X2 



+ + "r-. ^ + 



bx'i bx'i bx'i ax'i bx'i abx'i bx'i 
2?/^X2 ^ y^ X2 ^ XX2 ^xx2 yxx2 yxx2 ^ 



abx'i abx'i h^x'i bx'i IPx'i bx'i 
2y'^xx2 yx'^X2 yx^X2 x^ X2 yx2^ y"^ X2^ 



+ 



abx'i b'^x'i abx'i h'^x'i abx'i abx'i 
XX2' ^yxx2' x^ X2^ 8x'i ^ 4xx'i ^ 4x2x'i 4x2x'i 



6^ x'l abx'i IP' x'l a^ a^ a? a? x 

2xX2x'i 8x'2 4xx'2 4X2X'2 4X2X'2 2XX2X'2 



+ — ^ + 



dp aP aP aP aP x aP 

2 y a;'2 2 X x'2 2 X x'2 2 X x'2 yx2x'2 y'^X2x'2 



ax'i ax'i bx'i abx'i abx'i abx'i 

XX2X'2 XX2X'2 2 ?/ X X2 x'2 X^ X2 x'2 



6^ x'l abx'i abx'i abx'i 



(C5) 



{I - X2 - y){l + {I - xf) (x2 + yf + y^ 

~ 2? — \ \ T^; \ ' ^^^) 

x^(x + X2 + y — 1) X2 + y 

here a = 2f (1 - C12) - x - X2, b = 2f (1 - C12) + x'^ + x'2, x'l = ^{1 - Ci), x'^ = 

^(1 - C2), C12 = C1C2 + S'iS'2Cos02, where C12 = cos(A;i, ^2), 5'i,2 = sin(A;i,2, A = 
xy + XX2 (C2 - cos 02 j ■ 

Let us note that the combination ^-^^^-^^/^^^"'"'"(x, X2, + 2^ j is finite at X2 — > 

hmit. In this integral the value of X2 is fixed by delta-function in phase volume ()2ip : 
X2 = ^ (2 — 2(x + y) + xy{l — Ci)). Energy conservation law gives 0<X2<2 — x — ?/. 



APPENDIX D: VECTOR AND SCALAR 4-DIMENSIONAL LOOP INTEGRALS 

We introduce the following shorthand for impulse integrals (we imply real part in r.h.s.): 

_ , d% 1 _ / d'^ki k^i 



where we have used the short notation for the integral denominators 

(0) = kl - \\ 
(1) = (p - hf - M\ {!') ^{p-k- k,f - M\ 
(2) = (r - k^f - m^ (2') = (r + A; - k^f - 

The integrals with two denominators 

y 



Ji'2 — La, Joi — La + 1, Ji2 — L\ + 1 + 
Ji2' — La, Jii' — L\ — 1, Joi' — La + I + 



i-y 

X 



l-x 



Inx, 



J22' — I/A — 1 — Ljj, J02 — I/A + 1 — Lp, 



J02' = I/A + 1 - In(^), Jv2' = I/A + 1 - 



2-1/ 

1-y 



ln(2-y), 



where: La = In L/3 = In /? = m'^/M^. 

The integrals with tree denominators (we put M = 1 and introduce the notation A 



•I012 = 
Joi'2 — 

<Il22' = 

•I012' = - 



2|/ 



Le lnXl-ln^y + -L} + 2Li2 



1 — X 

-^2 + Inx In — h I/i2(2;) 



1-x 
1 

1^ 12 
1 



-\n^y + Li2(l - - ln/31ny |> , 



2(1-^) 



ln^2;, Joii' = - [-^^2(1 - a;) - ^2] , 



X 



J< 



022' 



2 /3 



11'2 



1-y 



-^^2(1-2/), 



where Le = In ^, ^2 = ^^2(1) = 7r^/6. 
We also need two integrals with four denominators: 

1 



0122' — 



1,2 ^ - ^ 



+ 26 



Now we consider the vector integrals with tree denominators: 

Jq12 = «012 + /3oi2, J^V2 = (P - ^)'' «01'2 + /?o1'2, 

Jq^ii, = p'* aoii' + k^" /3oii' , Jo^i2' = «oi2' + (r + /c)'^ /3oi2' , 
J(;,2 = P'' "ii'2 + r'^ /3ii/2 + k^" cn'2, J'^22' = «i22' + P122' + k^" C122' 

■^022' ~ '"'^ Q;o22' + A;''/?022'. 



The coefficients 0;^^, Pijk and Cijk are tfie following: 
Oioi2 — — {J12 — Joi) , 

y 

If 2 \ 
P012 = — I J12 — J02 {J12 — Jqi) , 

y \ y ) 

ttoii' — — [Jqi ~ Joi' + xJoiii) , 

X 

If 2 \ 

/3oii' = — I Jiv — Jqv {Jqi — Jov + xJqu') 1 , 

X \ X J 

OiOl'2 — Z iJl'2 — JqI') , 

1 — X 

Poi'2 — Z (2Joi' ~ Jl'2 — J02 — ^Joi'2) , 

1 — X 

OiQW = -JZ (22; ( Ji2' - Jq2') - {1 + z) ( Ji2' - Joi + zJqi2')) , 

[i-zy 

/3oi2' = ~ _ \2 (2 ("^12' - Jqi + zJqi2') - (1 + ^) (J12' - J02')) , 
(i z) 

OC022' = — {Jq2' — J02 + zJq22') ■, 

z 

P022' — — {J22' — Jo2') , 

z 

«ii'2 = Jii'2 + {J12 — Ji'2) , 

1-2/ 

Pll'2 = Z {Jl'2 — J12) , 

Cll'2 = {JlV - 2Jl2 + Jl'2 - (2 - y) Jll'2) , 

1 -y 

q;i22' — (Ji2' — J12) , 

1 -y 

P122' — J122' + -; (J12 — Ji2') , 

1-2/ 

C122' = :; (2Ji2 — Ji2' — J22' — yJ\22') ■ 

1 - y 

The vector integrals with four denominators: 

■^0^122' = "0122' + l3oi22' + (r + kf C0122'. 
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The coefficients aijk, jSijk and Cijk are the following: 



1 

Q;oi22' 



2(1 -x) 



^122' + -7-— {zJ022' - (1 - 2?/ + z)Jqi2> - y (Joi2 - ^^0122')) 

1-2/ 



/?oi22' = 7T71 TTl ^ [{z + y- zy- 1) Ji22' - z{l -2y + z)Jq22' + (1 - 2)^-^012' + 

2(1 — a;)(l — y)z 

+ (y(l + ;2) - 2^) (^Joi22' - J012)] , (DIO) 

C0122' = ^i(I^[(l/-2.)Jl22' + 

+ (-^"012 - Zjm22') - yzJo22' - (y(l + Z) - 2z) J012') 

1 - ?/ 
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y / X 
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0.8 


0.9 


0.9 


41.000 


8.278 


2.833 


1.250 


0.644 


0.371 


0.232 


0.156 


0.8 




33.111 


8.500 


3.333 


1.611 
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0.356 


0.7 






25.500 
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20.000 
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16.111 
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1.558 


0.4 












13.347 


4.875 


2.494 


0.3 














11.375 


4.364 


0.2 
















9.975 



TABLE I: The value of ^>}^ (i.e. Born inner bremsstrahlung part, see HA5() ). 



y / X 
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0.8 
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-0.020 


0.6 








-0.772 


-0.203 


-0.084 


-0.043 


-0.025 


0.5 










-0.521 


-0.140 


-0.059 


-0.031 


0.4 












-0.348 


-0.093 


-0.040 


0.3 














-0.221 


-0.056 


0.2 
















-0.118 



TABLE II: The value of ^ (L^ - 1) (see (jBTJ). 
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-2.362 


-1.978 


-1.596 


-1.229 


-0.994 


-1.161 


0.7 






-2.657 


-2.493 


-2.248 


-2.012 


-1.850 


-1.941 


0.6 








-2.716 


-2.600 


-2.438 


-2.333 


-2.431 


0.5 










-2.892 


-2.748 


-2.661 


-2.768 


0.4 












-3.137 


-2.979 


-3.049 


0.3 














-3.575 


-3.405 


0.2 
















-4.324 



TABLE III: The value of K{x,y) (see (jSSl). 
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